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Abstract 

For free fermions at finite density, the Pauli exclusion principle is responsible for 
the existence of a Fermi surface and the consequent presence of low energy spectral 
weight over a finite range of momenta. We investigate the extent to which this effect 
occurs in strongly interacting quantum matter with a holographic dual. We obtain 
the low energy current-current spectral weight in two holographic frameworks at finite 
density: systems exhibiting semi-local quantum criticality (with a low temperature 
entropy density vanishing like s ~ T n ), and a probe D3/D5 system. For the semi-local 
theory with < n < 2 we find a sharp discontinuity in the transverse spectral weight at 
a nonzero momentum k±. The case n = 1 is found to have additional symmetries and 
is soluble even at nonzero temperature. We show that this case exhibits a robust linear 
in temperature resistivity in the presence of random charged impurities. For the probe 
D3/D5 system we find an analytic expression for the low energy spectral weight as a 
function of momentum. The spectral weight is supported below a specific momentum 
k± and is exponentially suppressed at higher momenta. 



1 Introduction 



The Fermi liquid is a remarkable state of matter. A modern understanding of Fermi liquid 
theory is as a self-consistently robust organization of the low energy degrees of freedom 
in metals [TJ [2]- The scaling towards a surface in momentum space can be grounded to 
some extent in the Luttinger theorem that relates the existence of Fermi surfaces to the 
fact that there is a charge density [21 HIE]. ^ n a simple qualitative picture of a Fermi liquid, 
the exclusion principle plays a double role. By preventing the single particle degrees of 
freedom from condensing in their ground state, it both prevents spontaneous breaking of 
electromagnetic symmetry (at least at the natural scale set by the chemical potential) and 
furthermore pushes some low energy excitations of the system out to a finite momentum. 

In attempting to move beyond Fermi liquid theory, a crucial question is the extent to 
which the two roles of the exclusion principle can be decoupled. For instance, perhaps 
strong interactions persisting to the lowest energy scales can allow the charge-carrying low 
energy degrees of freedom to spread out over a volume rather than a surface in momentum 
space, or to collapse to a point in momentum space without triggering Bose condensation. 
In this paper we will explore this question in two different holographic setups. 

The momentum space structure of a Fermi liquid underlies much of its characteristic 
phenomenology. One should be careful, however. It is sometimes stated that due to having 
a surface in momentum space of low energy degrees of freedom, a Fermi liquid has infinitely 
many more degrees of freedom than bosonic systems. The Fermi liquid low temperature 
specific heat c ~ T is indeed larger than the c ~ T d specific heat of free relativistic bosons. 
However, more generally bosons have c ~ T d / Z , and so there is a quantitative sense in 
which bosons with dynamical critical exponent z = d have the same number of low energy 
excitations as a Fermi liquid. Similarly, the logarithmic violation of the area law in the 
entanglement entropy of Fermi liquids j6j [7] is closely tied to the decoupling of low energy 
degrees of freedom at different points on the Fermi surface [8] [9]. However, at least in 
holographic settings as well as Fermi liquids, this logarithmic violation seems to be tied to 
the system having a specific hyperscaling violation exponent, 6 = d — 1 |10l 111) . This is 
a kinematic feature of the system that a priori is logically independent of the momentum 
space structure of the degrees of freedom [12] . 

Other phenomenological features of Fermi liquids, in contrast, are tied unambiguously 
to their particular momentum space structure. Quantum oscillations in a magnetic field 
are a direct consequence of non-analyticity in momentum space [13] . The T 2 dependence 
of the resistivity of a clean Fermi liquid due to umklapp scattering depends crucially on 
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the interplay between the lattice and Fermi momenta, e.g. |14j . Friedel oscillations induced 
by an external charge depend on the fact that the static susceptibility of free fermions is 
non-analytic at momentum equal to 2k F , e.g. |15j . 

This paper will focus on the observables that directly characterize the existence of 
charged low energy excitations as a function of momentum. These are the spectral den- 
sities - imaginary parts of the retarded Green's functions - of currents and charges. The 
linear response of a medium to an electromagnetic field may be divided into transverse and 
longitudinal modes. The quantities of interest are 

We will compute the spectral densities^ above in two holographic theories, introduced below. 
Let us first recall the results obtained for a free nonrelativistic fermion in 2+1 dimensions. 
Throughout this paper we will consider 2+1 dimensional systems. The derivation of these 
formulae is sketched in appendix [A) 

a ± (k) 



jfe- k<2kF , *#) 

k>2k F 



k < 2k F 



kyJk 2 F -k 2 /4, . ( 2 ) 

k>2k F 

In both channels we see that there are no low energy excitations of the system above twice 
the Fermi momentum k > 2k F . In the longitudinal channel there is a strong discontinuity 
at k = 2kp, while in the transverse channel it is the derivative of the spectral weight that 
is discontinuous. The weaker discontinuity in the transverse channel has a geometrical 
origin. At low energies the fermions that contribute to the spectral weight are on the Fermi 
surface. If we now look at external momentum 2k F , the only low energy fermions that could 
contribute to the conductivity have momentum parallel to the external k = 2k F momentum 
and so do not contribute to the transverse conductivity. The main objective of this paper 
is to establish the extent to which a momentum space structure analogous to Q exists in 
known holographic theories. 

Before moving on, we would like to contrast the spectral weights with other physically 
important quantities, the static susceptibilities. These are defined as 

X ±(k) = hmG* X|| (*) = lnnGK^A;). (3) 

The static susceptibilities determine the screening of static currents and charges by the 
medium. For our free nonrelativistic fermion these can be easily obtained by performing 



Throughout this paper we will refer to ImG fl (w)/ ui as the spectral density. This quantity is always 
nonvanishing as ui — > at finite temperature, and arises naturally in many contexts. 
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the integrals in appendix [Aj The susceptibilities exhibit a kink at k = 2kp, but do not 
vanish for k > 2/cp. Upon Fourier transforming the screened Maxwell field propagator, 
this kink leads to spatial Friedel oscillations in the response to a point charge. While in 
a Fermi liquid the static susceptibility is dominated by the low energy degrees of freedom 
and exhibits a (weaker) non-analyticity at the same momentum as the spectral weight, in 
general the static susceptibility is not a direct probe of low energy degrees of freedom. In 
fact, from the Kramers-Kronig relations, the static susceptibility is equal to the spectral 
weight integrated over all energies. We will see explicitly in our first holographic example 
how the static susceptibility cannot be determined from purely low energy data. The focus 
in the remainder of this paper shall be on the distribution of low energy spectral weight in 
momentum space. 

Holographic descriptions of (compressible) quantum matter at zero temperature exhibit 
an emergent scaling symmetry in the far interior of the dual spacetime [TS]- This scaling 
symmetry is characterized by a dynamical critical exponent z [17] as well as a hyperscaling 
violation exponent 6 [T%1 [TB] . Because the scaling in all known holographic geometries 
is towards the origin in momentum space, u) ~ k z , the momentum space structure of the 
low energy charged degrees of freedom is determined by z. In particular, if z < oo is finite, 
then the spectral weig hts at any nonzero k are exponentially suppressed as ui — > |12j . 
In all of these systems with z < oo, a strongly coupled version of the exclusion principle 
does not appear to be active j^] And yet, these geometries can in general be stable against 
symmetry-breaking instabilities. These non-Pauli excluding geometries include those with 
logarithmic entanglement entropy jlOl 111] , 

In this paper we will obtain the spectral densities in two holographic frameworks 

that do exhibit a strong coupling remnant of Pauli exclusion physics. The first concerns 

theories with z = oo. In these systems, time scales at the IR fixed point, but space does not. 

It is therefore possible to have significant low energy spectral weight at nonzero momentum 

[HI [201 [12] . Such theories have rightly been termed semi-local quantum liquids [21, 22J; they 

are semi-local because while momenta do not scale, the scaling of correlators with energy 

is generally momentum-dependent. We will see examples of this in section [2] below. In a 

truly locally critical theory, the IR scaling would be completely independent of momentum. 

One of the main results in this paper will be the observation that momentum-dependent 

2 The discussion in this paragraph and throughout this paper refers to leading order in the bulk semiclas- 
sical expansion, corresponding to leading order at large N in the boundary field theory. This is the limit in 
which holography provides the greatest theoretical and conceptual simplification. 
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operator dimensions can lead to important features in the low energy spectral weight. Semi- 
local quantum criticality can also lead to strong dissipation in fermionic correlators [23j . 
Semi-locally critical IR geometries are parameterized by the quantity rj = —6/z which is 
kept finite and positive as z — > oo. They are conformal to AdS2 xl 2 . At low temperatures, 
they have an entropy density scaling like s ~ T' n . For r/ > the entropy density vanishes 
at zero temperature, avoiding the 'entropy catastrophe' of pure AdS2 x IR 2 . On the other 
hand, geometries with ij > have divergent curvatures in the far IR, possibly requiring 
additional physics to resolve. 

The second set of theories we consider are probe brane constructions at finite density. 
In these cases the charged degrees of freedom are parametrically diluted by neutral modes. 
In the bulk this means that the brane degrees of freedom do not backreact on the metric 
and therefore we do not obtain scaling geometries of the form discussed in the previous 
paragraph. This important physics thrown out in the probe limit is partially compensated 
for by the fact that the brane action is nonlinear in the bulk Maxwell field. At finite 
charge density, due to these nonlinearities, in the far IR the equations describing charge 
and current fluctuations become independent of momentum. This observation suggests that 
low energy spectral weight will exist at nonzero momenta, reminiscent of Fermi surfaces or 
semi-local criticality |24j . Unlike the scaling geometries of the previous paragraph, however, 
the IR does not decouple in a simple way and the full spacetime is needed to compute the 
momentum dependence of the spectral densities ([!]). This computation will be the topic of 
section El 

The Pauli exclusion principle is a statement about a quantum mechanical wavefunction 
written in terms of fermionic degrees of freedom. We are interested instead in strongly 
interacting effective low energy descriptions of finite density systems that do not admit a 
quasiparticle description, fermionic or otherwise. Indeed, it would seem to be the case in 
holographic frameworks that it is not a well defined question, from a universal IR point 
to view, to ask if the charge is carried by bosons or fermions. In particular, we should 
resist, perhaps, contaminating our IR perspective with information such as whether the UV 
charged degrees of freedom in an intersecting brane model are purely fermionic or not (e.g. 
D3/D5 versus D2/D6 etc.). In this paper we will find that strongly interacting 'IR soup' 
can nonetheless exhibit a feature characteristic of Pauli exclusion, namely the presence of 
low energy degrees of freedom at nonzero momenta. One is led to wonder about the extent 
to which Pauli exclusion may be a special case of more profound and general quantum 
dynamics. 
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2 Semi-local quantum liquids 



The ubiquity of semi-local quantum criticality in holographic theories at finite density is one 
of the more interesting facts to emerge from a condensed matter perspective on holography. 
While the simplest A1IS2 x M 2 geometry may be pathological due to having a macroscopic 
ground state entropy density, this is not the case for generic semi-locally critical geometries. 
It was emphasized in [12] that theories in which rj = — 9/z > is kept fixed as z — > 00 
have semi-locally critical spectral weights, with an entropy density vanishing like at low 
temperatures. These spacetimes are conformal to AdS2 x and appeared in previous 
works including |25|, [26l [271 128] . We will describe the embedding of semi-local criticality 



into string theory in a particularly interesting case (77 = 1) in section 2.3 below. 

Following [18] . using the notation of |llU12j . our starting point is the Einstein-Maxwell- 
dilaton theory in 3+1 bulk dimensions 

C = -^i? - -^Z($)F 2 - \{d<$>) 2 - 7TTr2 V (®) ■ (4) 
2k 2 4e 2 k 2k 2 L 2 

The low energy spectral weight at low temperatures can be obtained from the IR near- 
horizon geometry of near-extremal black holes in these theories. We therefore only need the 
forms of the potential and gauge- kinetic function to leading order as $ — > 00, as the dilaton 
will diverge in the far IR of these solutions. These asymptotic behaviors are taken to be 

Z(*) = r 2 Z 2 e^ , V{*) = -r 2 V 2 e-^ . (5) 

The factors of Tq will be convenient shortly. As presented above, Zq and Vq have units of 
inverse length. They will carry the units of momentum in the semi-locally critical theory. 
This theory admits a scaling solution where the parameters a and (3 will be related to the 
dynamical critical exponent z and hyperscaling violation exponent 9. Our interest is in the 
limit z — > 00 and — > —00 with 77 = — 9/z held fixed. This will be seen to correspond to 

a = -P= . 2 • (6) 

Black hole solutions to these theories may be found explicitly and take the form 



ds 2 = (-f(r)dt 2 + g(r)dr 2 + dx 2 + dy 2 ) , A = —h(r)dt, * = 4/1 + -log— . (7) 
r z ' k V V r o 

Evaluated on the logarithmically running dilaton, the functions of the dilaton in the action 

become Z(Q) = Z 2 r 2 and V(&) = —V^r 2 . The functions /, g and h are given by 



where the constants are 



V 2 



90=772 1 + " . h ° = v fTTL — ' ■ ( 9 ) 
and the emblackening factor \ is 

/ N 2(1+1/*)) 

X(r) = !-(-) ■ ( 10 ) 

Here is the horizon radius, which determines the temperature of these black hole solutions 
to be 

T=(l+ 1 -)^ = %-^. (11) 



From these expressions, we see that the entropy density s indeed vanishes at low tempera- 
tures: s oc r^ 2 ~ XT At zero temperature, x = 1, the metric transforms covariantly under 
the scaling t — > Xt, r — Y X v ^ 2 r, ds — > A _J " 2 <is. This is our local criticality, time scales but 
space does not, together with hyperscaling violation parameterized by rj. 

To obtain the current-current correlations functions, we need to perturb the solutions 
above. Taking the momentum of the perturbations to be in the x direction, the perturba- 
tions will decouple into modes that are even or odd under y f-> —y. Thus upon perturbing, 
and picking the gauge A r = g^ r = 0, we will obtain coupled equations for the transverse 
modes {5A y , 6g xy , Sg yt } and the longitudinal modes {6A t , SA X , 5g tt , Sg xt , Sg xx , 5g yy , 5$}. The 
perturbations are taken to have the form of a function of r times e - lujt + tkx . 

Due to a residual gauge symmetry, the equations for the perturbations above will be a 
mix of first and second order differential equations. A set of independent and purely second 
order differential equations are obtained by introducing gauge-invariant variables. These 
variables are derived in appendix [B] and can be taken to have the following form. In the 
transverse channel 

k r 2 
Tp±,i = 77J^ A y ' V>±,2 = j2^ uS 9xy + kSgyt) , (12) 

and in the longitudinal channel 

^ = & SAt ~ Z Q L^v X (r) 5Stt > ^ 

^ = V 5 ** + 2^ (1 + „) (1 + 2 P) - (1 + 2 V ) X (r) dr [ ^ ^ + 5g ™ )] ' (W) 
5<& r 2 
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We have defined k = \J\ + l/rjk/Vo. From the perturbed Einstein-Maxwell-dilaton equa- 
tions of motion we obtain coupled equations for each set of the gauge invariant variables. 
These equations can be further decoupled to a large extent. In the longitudinal case espe- 
cially, the equations are somewhat messy - the details can be found in appendix [C| 

2.1 Transverse channel 

The transverse channel at zero temperature was decoupled in [T2]. Here we extend the 
decoupling to nonzero temperatures. Taking our cue from [TH], we define 

9?± ' ±(r) = r 4/^2 _ k 2 x{r) [P±(rW±,i(r) + ^L )2 (r)J , (16) 

where we have rescaled if)±i = Zq^JI + Tjij)± i and the functions p± are given in appendix 
[Cj In terms of these new variables, the equations of motion are decoupled 



= Pt,±(r)+ l A% K i V 2 ( \ + ^ ' '\ r 2 2 2 21 f ^' /J 1 ^,±0)' ( 17 ) 
where q± (x) is given by 

ff±(x) = 4(1 + V) 2 ~ 8(1 + v) (3 - V ± 2^1 + 2^ x - 15r/ 2 x 2 • (18) 

At nonzero temperatures we are only able to solve these equations in the case r] = 1 (as 
well as the AdS2 x R 2 case r/ = 0, solved in |14j). to be discussed in more detail below. 



( x (r)) , 4(1 + T?) 2 [r 4 /^ 2 - k 2 X (r)] 



At T = we can quickly recover the results of [12]. At zero temperature x = 1 i n (17 



18). The resulting equation is solved in terms of Bessel functions for all r\. The solution 



satisfying infalling boundary conditions at the horizon (r — > 00) is 



vlAr) = , (19) 



where the indices are 
-i- 1 



k\ 2 ., , . N / A; x 2 



^I = 2M 5 + 7 ?( 2 + 7 ?) + 4 ( 1 + r ?) 2 ^J ±4(1 + 77)^1 + 2(1 + 77)^-] . (20) 

The retarded Green's function is now found by expanding the solution near the boundary 
r — > |29} I12j. From the asymptotic expansion of Hjy(x) for x->0we obtain 

V±,±(r) oc y^(r- 2 ^ + <?f,±(u;,£;)r 2 ^) , (21) 
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with 

Sf,±(w,fc) oc w 2 ^W. (22) 

This scaling pertains to both the real and imaginary parts of the Green's function. 

We are now in a position to obtain an interesting result for the transverse spectral 
weight. We have obtained the Green's function in the IR near- horizon geometry, and have 
not specified how this is connected onto an asymptotically AdS± or other spacetime which 
would define the UV of the system. Recall, however, that to leading order at low frequencies, 
the imaginary part of the full UV Green's function is given by a sum over the imaginary 
parts of the IR Green's functions of all operators that overlap with the UV operator [30]. In 
the present context, the exponent v7 leads to the strongest IR behavior and so we obtain 
that at small frequencies 

lmG R ±Jx (oo,k) oc Imgf_(ui,k) oc w 2 "!^' . (23) 

And therefore the transverse spectral weight ([I]) is found to be 

ImG 1 } T (u,k) oo k < k+ 

a ± (k) = lim J±J±V ; = I . (24) 

w I k > K 

The critical momentum is given by 

g = ^ g(W|^p EES, ( o<„< 2) . (25) 

Outside the range < rj < 2, the transverse spectral weight is zero everywhere. 



The reader that finds (24) unpalatable, or worries that dividing by some given power 
of uj is artificial, might like to turn on a small temperature. While we have not solved 
the equations in general at finite temperature, it is straightforward to show that at low 
temperatures and for uj <C T one has 

ReQ R (uj,k) ~T 2u , lmg R {uj,k)~ujT 2v - 1 . (26) 

This expression is true for all r\ in both the transverse and longitudinal channels. It follows 
in particular that 

a ± (k) oc T 2 ^- 1 . (27) 

This expression is finite and nonzero for all momenta. As we take the low temperature 
limit, the spectral weight diverges for k < k± and goes to zero for k > k+. 
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The result (24) for the zero temperature transverse spectral weight, or the thermally 



regularized version (27), dramatically reveals the existence of low energy spectral weight 
over a nonzero but finite range of momenta. In our theories with < r/ < 2 - this does not 
include the simplest case of AdS2 x M 2 - a strongly interacting cousin of Pauli exclusion 
appears to be operating. The physics is tied in a precise way to the emergence of semi-local 
criticality [21] . For k < k* the leading IR operator excited by the transverse current is 
relevant under the local scaling, while for k > k+ it is irrelevant. This follows from the fact 
that the dimension of the IR operator is A = v + |, and so becomes relevant if v < I. 
But v < \ corresponds to the range of momenta over which the zero temperature spectral 



weight diverges in (24). 



One must be careful not to equate the k = values of (22 ) and (27) with the optical and 
d.c. conductivities. For the conductivities, it is crucial that the limit k — > be taken before 
the uj — > limit. At exactly k = 0, the two gauge-invariant transverse modes decouple. 
The transverse current operator at k = is found to correspond to the index i/i"(0), and 



not z^_|_(0) as one would obtain by naively evaluating (22) and (27) at k = 0. The optical 
conductivity is therefore while the critical contribution to the d.c. 

conductivity is a± c . ~ T 2+7? . These expressions agree with previous results in e.g. [26]. The 
critical d.c. conductivity is swamped by a translation-invariance delta function. In section 



2.3 below we discuss the effects of impurities that render the d.c. conductivity finite. There 



we will obtain a strong d.c. resistivity, that is linear in temperature when 77 = 1. 



The critical momentum fc* is set in (25 ) by the data {77, Vb} that appear in the definition 
of the bulk theory and the IR scaling solution. It is natural to consider the dimensionless 
ratio of the critical momentum divided by the (square root of the) total charge density 
of the system. The total charge density is defined in the UV but, in the absence of bulk 
charges, can also be obtained from the IR because Gauss's law implies that the integrated 
flux is independent of radius [31]. The charge density is then 

„ = /,(Z(,)^__£^_ W* (28) 

The dimensionless ratio is therefore 

kl = Vp V (2- V ) + 2^(1 + y)(2- 77) 

Q Z 4(l + r?;)3/2 • 1 ' 

This ratio can be obtained from the parameters characterizing the IR Lagrangian, i.e. it is 
a property of the IR theory, but does not appear to be a universal number. 



We have found the low energy spectral weight ( 24 ) using only IR data about the theory. 



As we anticipated in the introduction, this is not the case for the static susceptibilities 



Recall that, schematically, the full zero temperature Green's function at low frequencies 
takes the form 



G 



R 

J±J 



(U3,k) 



A(k) + B(k)G^_(u,k) 



(30) 



c(k) + D(k)g*_(oj,k) ' 

Here {A, B, C, D} are real functions of momentum that are determined by solving the bulk 
differential equations in the 'far' region, away from the horizon, at u) = 0. The full structure 
is in fact more complicated than this, due to mixing among different modes between the IR 
and the UV [30J, but the above expression is sufficient for illustrative purposes. Because 



is ± (k) > in the IR Green's function (22), when we take lo — > in the full Green's function 



to obtain the susceptibility, the contribution from the IR Green's function will go to zero. 
We are left with x±(k) ~ A{k)/C{k). This function of momentum is not determined by 
near-horizon IR data, but requires solving the Einstein-Maxwell-dilaton equations in the 



full geometry. It follows that ( 24 ) cannot be directly related to the existence of (transverse) 
Friedel oscillations. 



2.2 Longitudinal channel 

The equations of motion describing the three gauge-invariant longitudinal perturbations 
(13) - (15) are given in appendix [Cj It is possible to obtain decoupled equations for two 
linear combinations of the three modes, but the resulting equations are still somewhat 
messy. Here we present results at zero temperature, while below we will solve the equations 
at finite temperature in the case rj = 1. 

For T = 0, the longitudinal equations given in appendix [C] simplify to 







<P\\ 



4(1 + T7) 2 
r 2 r] 2 VQ 



fc2 )^||, 



(31) 



Here ipu^ (i = 1,2,3) are linear combinations of the gauge-invariant modes (13) - (15) and 
are defined in the appendix. The matrix M has constant entries given by 



M 



/_ 3 4_ 

4 2+r) r) 
8(l+r ? )(fc 2 -l) 

»7(2(l+r/)fc 2 -ri 

4(1+??) 
\ ??(2+r?) 



9+14?? 

,2 





(2+r?)(2+3r?) 
4?? 2 





16(l + T?) 2 fc 2 \ 
n 



16(l+??)fc 2 
r?(2(l+r?)fc 2 -r?) 
(2+r?)(2+3??) 
4r? 2 



(32) 



By diagonalizing M we obtain three decoupled differential equations. These can be solved 
in terms of Bessel functions. Imposing infalling boundary conditions at the horizon r — > oo 
gives the solutions again in terms of a Hankel function 

V j " 



(i) 



V||,a(r) = VrHL 



(33) 
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Here a = {0, +, — } labels the decoupled modes and the indices are given by 

^(k)= 1 -V L \k + 4l)\ (34) 



Similarly to the transverse channel, from the asymptotic expansion of the Hankel function 
we find that the Green's functions have scaling 

Gf a (u,k) oc u? v \\ . (36) 



Unlike in the transverse case, we see in equations (34) and (35) that for all r\ and 
momentum, the indices Uvf > |. It follows that the longitudinal low energy spectral weight 
([I]) always vanishes. It does not exhibit the dramatic divergence below a critical momentum 



that we found in (24) for the transverse sector. Relatedly, the finite temperature low energy 
spectral weight (27) will now be uniformly small at low temperatures over all momenta. 
In terms of semi-local criticality, this is the statement that the longitudinal operators are 
irrelevant in the IR theory for all momenta. 

It is interesting that these systems, with < ij < 2, exhibit a sharper structure in 
momentum space in the transverse channel than in the longitudinal channel. This is in 
contrast to the free fermions that we reviewed in the introduction. It is perhaps natural 
that the longitudinal channel operators are irrelevant in the semi-locally critical picture. 
The charge density at k = cannot be relevant: that would imply we had not isolated the 
correct finite density IR fixed point. It seems reasonable, in addition, to expect that the 
scaling dimension will increase with k. The irrelevance of the longitudinal channel operators 
for all momenta further implies that our semi-locally critical IR theories are robust under 
addition of an ionic lattice El . 



2.3 The magical r\ — \ case: enhanced symmetries and linear resistivity 

In this section we will find that the case with rj = 1 has some interesting mathematical and 
phenomenological properties. This is perhaps made more interesting by the fact that this 
case arises easily from string theory compactificationsj^] Extremal black holes in AdS± x S 7 
with one of the four U(l) charges being zero all have a near horizon geometry that is 
locally critical with r\ = 1 [32]. These critical geometries therefore describe finite density 



We thank Aristomenis Donos for pointing out to us that these embeddings all realize r) — 1. 
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phases of the world volume theory of N coincident M2 branes [33]. Furthermore, if the three 
nonvanishing charges are equal, the background solutions can be obtained from a truncated 
theory with one gauge field and one scalar, of the type we are considering in this paper (e.g. 
|25j , note however that finite momentum perturbations will probably require a larger set of 
fields due to an F A F term in the action that sources axionic fields |32j). 



For 77 = 1, the transverse equations (17) reduce to 



„ ( 16 (r^ 2 -k 2 X (r)) 16-15 X 2 (r) 8(l±vT + g 
° = ^< ±(r) + [ V^xHr) + 4rVM rW) ' * ±M) ' ( ' 

This value of r\ is the only nonzero value for which we have been able to solve the finite 
temperature equations analytically. This equation can be solved in terms of hypergeometric 
functions. The solution satisfying infalling boundary conditions at the horizon is given by 

¥>_L,±(r) = r^xir^-^r (T(c)T(d)T(l - uf) (r/r h ) 2u i 2 F X (c, d, 1 + v^/ri) 

+ K -»f}) , (38) 



where c = —~^f + d = 1 — + and the indices are those quoted in section 
setting 7] = 1: 



2.1 



f = V2Vl + k 2 ± \Jl + 2k 2 , (39) 



and we used the temperature T = Vb/(47rr^), from (11) with rj = 1. 

Expanding near the boundary, as in previous sections, we obtain the retarded Green's 
functions 

GZ^k) oc ^ J " T ^ . (40) 

1 ^ 4ttT 2 / \ 47rT 2 J 

Expanding this expression at small frequencies gives 

ReQf ± {u,k)~T 2u i , Imgf ± (u,k) ^ uT 2 ^ 1 . (41) 



Here we have recovered explicitly the general result ( 26 ) that holds for all values of 77. 



The Green's function (40) is of a form that strongly suggests the existence of hidden 
SL(2, R) symmetries. More precisely, it is the form that follows from an SL(2, R) x SL(2, R) 
symmetric system in which a hyperbolic generator has been 'level matched' |34| . It seems 
likely that this symmetry is directly related to the AdS% factor that appears in the string 
theory uplift of these geometries [25J. See [35] for a related earlier appearance of AdS^ 
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in these types of solutions. This offers the hope of obtaining an explicit dual to our IR 
solutions and thereby a microscopic understanding of the hyperscaling violation and the 
nature of the momentum dependence of the spectral weight. 

Turning now to the longitudinal channel, at finite temperatures and with rj = 1, the 
equations for tp^ in appendix |C| simplify to 



M, 



+ 



Ni. 



r 2 x(r) r 



r 2 x 2 (r) 



16 (rW - k 2 X {r)) \ 



(42) 



where the matrices M and N are given by 



M 



( 



\ 



2 \ 



-f -64A: 

4(2fc 2 -3) o 16(l-2fc 2 ) 
2fc 2 +l l+2fc 2 

-I -8 



/ 1 

4 



N 











V 



1+2A; 2 









63 16(3fc 2 -l) 
4 2fc 2 +l 
1 
4 



(43) 



Magically, it seems to us, these two matrices commute so they can be simultaneously diag- 
onalized. The diagonalized equations can then be solved analytically in terms of hypergeo- 
metric functions, rather similarly to the transverse case. The infalling solution is 

<P\\Ar) = r^xCr)™ (r(c)r(d)r(l - u^) {r/r h f v \\ 2 F 1 (c,d, 1 + ^,r 4 /rf) 

+ {v§ -> . (44) 



Here c 



and d 



-> a - n 2~ ail,a " _ _s a ~ 5fT 2^' As a bove, a = {0, +, -} labels the 

decoupled modes, and sq = 1 while s± = 0. The indices are obtained by setting rj = 1 in 



our previous expressions ( 34 ) and ( 35 ) 



l + 2k 2 
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f 2k 2 ± - A /l + -k 2 
3 3 V 2 



(45) 



Once again imposing ingoing boundary conditions at the horizon and expanding near 
the asymptotic boundary, we obtain the retarded Green's functions 



g?Ju,,k)<xT 



2i/f? ' 



'AnT 



+ 



1+1/? 



+ s a r 



^aJ I I 
"4ttT "t 2 



r 



"4ttT 



+ 



+ s r 



to) 
'4vrT 



l-i/? 



(46) 



These are again found to have the form expected from an underlying level matched SL(2, M.) x 
SL(2,M) symmetry. The transverse Green's function (40) is also in the form of the above 



equation (46) with s± = g. Comparing these Green's functions with those derived from an 
SL(2, M) x SL(2,M) symmetry in |34| . we see that the modes live in a symmetric lowest 
weight representation of 5L(2,M) x SL(2,M) with weight A = (v + l)/2 and the matched 
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hyperbolic (i.e. H — K generator in the SL(2, M.) algebra) momentum is, in the notation of 
HH, x = 2s. 

In possession of the finite temperature longitudinal spectral weight, we can obtain an 
interesting phenomenological result. In particular we consider the d.c. electrical resistivity. 
To obtain a nonzero resistivity it is crucial to break translation invariance. One way to do 
this is via random impurities. For a locally critical theory in the presence of random charged 
impurities, the density-density spectral weight ([!]) controls the d.c. electrical resistivity 
through the formula [T11I361I37] 

d 2 k l2l . lmGj tJt (u,k) 



I 



Um ImG„Kfc) ^ f dkk s T 2^i (47) 

w^O CO J 



(27T) 

The function v7 {k) is monotonically increasing with k. At low temperatures the integral is 
therefore dominated by small momenta. Up to corrections that are logarithmic in temper- 
ature, the leading low temperature behavior of the resistivity is found to be 

r imp . -T 2 ^- 1 = T". (48) 

Curiously, we see that precisely the value rj = 1, which we have just seen to have simple 
string theory realizations and enhanced symmetries relative to the other cases, gives a linear 
in temperature resistivity. A linear in temperature resistivity is widely observed in quantum 
critical metals, see e.g. [38] . 

The density-density spectral weight also controls the d.c. resistivity due to umklapp 
scattering by an ionic lattice. It was shown in [TJ] that in a locally critical theory, umklapp 
scattering leads to a resistivity 

r umk iap P ~ T 2 ^- 1 , (49) 

where ki is the lattice momentum. This result is UV sensitive, as the power itself depends 
on lattice data. In contrast, the impurity resistivity ( |48[ ) only depends on the IR theory. 
However, it is worth notin^Jthat the functions v7(k) are quite flat at small momenta, with 
the correction to the k = value going like k 4 . It follows that there will be a range of lattice 
momenta for which the d.c. resistivity due to umklapp scattering will also be very close to 
linear in temperature when 7] = 1. 

3 The D3/D5 system at finite density 

In this section we obtain the momentum space distribution of the low energy degrees of free- 
dom in a finite density holographic system with a known supersymmetric field theory dual. 



4 We thank Diego Hofman for emphasizing this point to us in discussions. 
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The field theory will be N = 4 SU (N) SYM theory at large N coupled to a single massless 
hypermultiplet in the fundamental representation. The hypermultiplet is furthermore con- 
strained to propagate on a 2+1 defect in the full 3+1 field theory dimensions. This system 
is known to be holographically dual to a single probe D5 brane wrapping AdSi x S 2 in the 
AdS& x S 5 background of type IIB string theory [391 HUl 01] . We will consider the system 
at finite density under the 'baryonic' global £7(1) symmetry that is dual to the Maxwell 
symmetry of the bulk probe D5 brane. There is a large literature considering Dp/Dq sys- 
tems in the probe limit at finite density, a slower exposition of the setup we consider may 
be found in [12]. Our results will hold for many other brane systems beyond the explicit 
supersymmetric D3/D5 theory we have chosen to study. 

The background spacetime is the Schwarzschild black hole in AdS§ x S 5 

ds 2 = L 2 f A {-f{r)dt 2 + AA + dx 2 + dy 2 + dz ^\ + d ^ _ (5Q) 

Here the emblackening factor /(r) = 1 — r 4 /r 4 . The temperature of the black hole is given 
by T = I/71T4-, see e.g. |29j . The probe D5 brane is described by the DBI action 

S D5 = -e-*T D5 J (fa^-det (*g + 2na'F) . (51) 

Here *g is the pullback of the metric to the worldvolume, F = dA is the field strength of 
the worldvolume Maxwell field. The prefactors are the string coupling and the -D5 brane 
tension. For the configurations we consider, the Wess-Zumino term in the brane action, 
and corresponding equations of motion, is zero. Taking the embedding AdS± x S 2 in the 
background spacetime, corresponding to massless fundamental fields, all that is necessary 
to place the system at finite density is to solve for the electrostatic potential At(r) on the 
worldvolume. The equations of motion are easily solved to give 

4(r) = --=£=. (52) 



'1 + Q 2 r 4 
:he 't Hoo: 

L 4 = Ac/ 2 . The constant of integration Q is proportional to the charge density of the 



Here we have introduced Q 2 = ^j-Q 2 , where the 't Hooft coupling A> 1 appears through 



-4 _ W2 



system. This can be seen by expanding (52) near the boundary at r = 0. More convenient 
for us will be the chemical potential 

In the probe limit, the excitations of the Maxwell field do not backreact on metric modes. 
This is technically a great simplification compared to the case we considered in the previous 
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section. Only the following perturbations are necessary: 8At,5A x ,5A y . Taking each mode 
to be a function of r times e - twt + lk v ^ Q ne finds the equation for the transverse mode 

where <5^4j_ = 5A X . The equation for the longitudinal mode 5 An = uj5A y + kSAt is 
6g 2 A; 2 r 3 / 2 -6J 2 (l + QV)(2Q 2 r 3 / + (l + 0V)/ , ) 1 



SA " ~ " ; :A '\+[ ^ - 7 — k -^— ) <U„ ■ - o 



/(1 + Q 2 r 4 ) (£; 2 /-u; 2 (l + Q 2 r 4 )) "V/ 2 / 1 + Q 



r 

(55) 

This decoupling of charge dynamics from heat flow and momentum conservation comes at 
a physical cost. For instance, the translation-invariance delta function in the conductivity 
is not visible in the probe limit |43j. 

We are primarily interested in zero temperature. This corresponds to setting / = 1 
in the above equations. In this limit we recover the equations in |42j . Our objective in 
this paper is to understand the momentum dependence of the current spectral density in 
these systems at vanishing frequency. We have managed to obtain this spectral weight 
analytically for the transverse modes. A low frequency expansion is implemented by setting 



SA± = e lujr [5AY_> + u6A£' +■■■ . (56) 

It is necessary to extract the infalling near-horizon behavior as we have done here to allow 
for the fact that the to — > and r — > oo limits do not commute. Said differently, the factor 
of e tU)T removes the irregular singular point in the equation at infinity. To zeroth order in 
uj we then find that the solution that is regular at the horizon is 



5A^\r) = sinh 



iarcsinh( (-lf^jQr) ,-l) + r ( 1 / 4 ) 21 



4x/7T 



(57) 



Here F is the incomplete elliptic integral of the first kind. The normalization of this solution 
is immaterial. The first order solution that is regular at the horizon can be written as 



\ _|_ Q2 r 4 r 

5Af{r) = c 5Af ] (r) - i ( r + ^ — ^ — ^ ) 5A { p (r) . ( 58 ) 

Here c is a constant. Obtaining this unexpectedly (to us at least) simple form of the 
first order solution required some nontrivial integration by parts of an integral expression 
that is obtained upon solving the differential equation for SA P . It may also be verified 



directly by plugging the form (58) into the differential equation for SA ^ and then using the 
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differential equation for 6A±. It is now straightforward to expand the full solution near 
the boundary r — > 0, giving 5A± = SA± n n + SA±^ n r + • • • . The retarded Green's function 
is given, as usual and up to an overall normalization, by the ratio of the normalizable and 
non-normalizable terms 



<^4±,n.n.(^,*0 ' 

The desired transverse spectral weight (dissipative conductivity) is then found to be 



lim 



lmGl Jx (u,k) 



a ± (0) ^csch— 



Here we introduced the zero temperature limit of the chemical potential (53) 

Q 4r(5/4) 2 



Mo 



lim a 



7T 



(60) 



(61) 



The exact expression (60) is the main result of this section and is plotted in the following 



figure [TJ Note the absence of hats on Q in the final expression ( 60 ) . We have verified this 
result by also computing the spectral weight numerically. The main physical point to take 




Figure 1: Spectral weight (dissipative conductivity) of the transverse current at zero fre- 
quency as a function of momentum for the D3/D5 system at finite charge density. This is 



a plot of equation ( 60 ) 



away from (60) is that for 



fc<fc* 



Q_ 

2/Uo 



(62) 



the system has spectral weight at low energies. For momenta greater than this, the low 
energy spectral weight is exponentially suppressed. 
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It has already been emphasized in [23] that gauge field fluctuations of the DBI action 
at finite density and zero temperature have the property that the differential equations 
describing the fluctuations become independent of momentum near the horizon. We can 



see this by taking r — > oo in (54) and (55), with / = 1 at zero temperature. It was explained 
that this implies the existence of low energy excitations at finite momentum. Our result 



( 60 ) makes this intuition explicit and shows furthermore that these excitations are present 
over a precise range of momentum < k < k*. 

The plot in figure [T] is satisfyingly along the lines of what one might expect for a strongly 
interacting quantum liquid. The well-defined vanishing of the free fermion transverse low 
energy spectral weig ht @ has been smoothened out by interactions. Nonetheless, the basic 
'exclusion principle' feature of low energy degrees of freedom extended in momentum space 
remains. 

Static screening of electrical currents is determined by the real rather than the imaginary 
part of the transverse current Green's function. The static Green's function is also com- 
putable analytically from our solutions above for 5A^ and SA^ . We obtain the transverse 
susceptibility 

XL(fc) = HmG« 4 (^) = -| X± (0)| ^coth^. (63) 

The magnetic response to a line of current is then obtained by Fourier transforming the 
screened boundary transverse gauge field propagator (k 2 — x±(k)) 1 • Curiously, using the 



susceptibility ( 63 ) we find that all the poles of the screened propagator are at purely imag- 
inary momenta. It follows that the screening leads to a purely exponential falloff of the 
magnetic field with distance from the current, with no oscillations. Possibly relatedly, no 
quantum oscillations in the magnetic susceptibility are observed in this and similar probe 
branes systems as a function of the magnetic field, e.g. [42j . 

For the longitudinal channel and at finite temperatures, we have not been able to solve 
for the spectral weight and static susceptibility analytically. In particular the zeroth order 
longitudinal equation with oj = does not appear solvable even at zero temperature. The 
relative complexity of the longitudinal equation is partially due to the fact that it supports 
a zero sound mode |441 14"2] . This feature dominates the spectral density at any small finite 
uj and k and is not a direct probe of structure at finite momentum. At zero temperature we 
solved the equations for 5A\\ numerically in a perturbative expansion in u, analogously to 



the expansion in ( 56 ) . This expansion pushes the zero sound feature down to arbitrarily low 
momenta. From the solution to the longitudinal equation we may obtain e.g. the retarded 
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Green's function for density fluctuations from 



k 2 5A ln Xu,k) 
k 2 — uj 2 5Aii D n (u, k) 



The longitudinal susceptibility 



X\\{k) = limG^ tJt (u,k), 



(64) 



(65) 



is found to have a qualitatively similar functional form to the transverse susceptibility (63), 



but with the opposite sign. We have verified numerically that, similarly to the transverse 
case, the screened propagator, which is now (k 2 + > nas poles at purely imaginary 

momenta. It follows that point charges are exponentially screened in this theory and do 
not cause Friedel-type oscillations. 

The spectral density is found to vanish with a higher power of uj than in the transverse 

(3) 

channel. We had to work numerically to order 8A\ to extract the leading nonvanishing 
spectral weight. Our numerical result for 

k 2 ImG* Jt (uj,k) 



3F„(fc) 



lim 



UJ 



(66) 



is shown in figure [2J Note the extra factor of k 2 /uj 2 in the definition of cry relative to 




Figure 2: Spectral weight (66) of the longitudinal current at zero frequency as a function of 



momentum for the D3/D5 system at finite charge density. This quantity has been obtained 
numerically. The zero sound feature has been pushed to arbitrarily small momentum. 



the density-density spectral weight defined in (JXJ) . That is, cry = cry (uj)k 2 /uj 2 . We see 
therefore that at a given finite momentum there is much less low energy spectral weight in 
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the longitudinal channel than in the transverse channel. At finite temperature, the fact that 
crii (w) ~ co 2 will translate into the spectral weight vanishing like T 2 at low temperature. 
The above said, the spectral weight that does exist is seen in figure [2] to be constrained to 
lie in a finite band of momenta. The dependence on momentum is qualitatively similar to 
the transverse channel shown in figure [I] and equation ( 60 ) . 

A quantity that can be extracted analytically at any temperature and for both channels 



are the momenta analogous to (62), above which the spectral weight is exponentially sup- 
pressed. As the temperature is increased, we expect to see spectral weight move to higher 
momenta as states are Boltzmann populated. The exponential decay at large momenta can 
be extracted with a WKB computation, cf. |45j . The imaginary part of the retarded Green's 
function is given by the probability for the excitation to tunnel from the boundary to the 
horizon. We can perform this computation directly at zero frequency, cj = 0. Mapping the 
equations to Schrodinger form and using k — > oo as the WKB parameter, one finds for both 
the longitudinal and transverse modes 

lmG R (k) ~ exp{ -2k T + dr } (67) 




In the low temperature limit, T — > 0, we recover the critical momentum in (62) separating 
regions with and without spectral weight. The same rate of exponential falloff with momen- 
tum is therefore seen to apply to both the transverse and longitudinal channels. At high 
temperatures the spectral weight is seen to be non- negligible for momenta (cf. [45j ) 



This is precisely the anticipated Boltzmann occupation of states. The hypergeometric 



function in the exponent of (68) curtails the decrease of k* with temperature so the finite 



region of momentum space in (62) remains active at zero temperature. 



4 Discussion 

Our starting point has been the observation that the existence of low energy degrees of 
freedom at nonzero momenta is a fundamental organizing principle for the best characterized 
compressible state of matter: the Fermi liquid. Given the central role this fact plays in the 
phenomenology and theory of metals, the extent to which such momentum space structure 



20 



appears in holographic phases of matter is currently an active topic of investigation 
fl~2"j RTF] . In this paper we have found the first holographic identification of a strong non- 
analyticity in the low energy spectral weight as a function of momentum that is visible 



to leading order in the bulk gravity limit (24). In a second class of holographic theories 
we have found low energy spectral weight that is exponentially localized below a critical 



momentum (60). 

Semi- local quantum criticality is among the most interesting new organizing principles to 
emerge from holographic studies of finite density matter. This criticality applies to IR fixed 
points in which time scales but space (and hence momentum) does not. The crucial extra 
ingredient of semi-local criticality is that the weights with which operators scale in time 
are momentum-dependent. Applied to current correlators, which are the quantities that 
universally detect the charge dynamics in the system at leading order in the bulk gravity 
limit, the momentum-dependence of the scaling dimension allows a sharp momentum-space 
structure to emerge at zero temperature. While the local scaling itself kinematically allows 
low energy spectral weight at all momenta, this spectral weight only participates in charge 
dynamics over a finite range of momenta for which the associated current operators are 
relevant in the IR. This effect has an interesting resemblance to the momentum space 
structure of the spectral weight of a Fermi liquid. 

Initial work on semi-local quantum criticality in finite density holography revolved 
around AdS 2 x R 2 [23]. It was emphasized in |I2J that this is a special case of holographic 
semi-local criticality, which is more generally parameterized by the quantity r/ discussed in 
the introduction. For rj > there is no ground state entropy density. For < r\ < 2 we have 
found an interesting momentum space structure. We have furthermore noted that rj = 1 is 
easily embedded into string theory and thereby realized as a finite density phase of the M2 
brane CFT dual to AdS A x S 7 . The same case of rj = 1 was found to exhibit mathematical 
simplifications and a linear in temperature resistivity in the presence of impurities. 

The key question facing semi-local criticality seems to be its generality. Can it arise 
naturally in non- holographic systems? Is it somehow an artifact of the large iV limit, 
perhaps describing multiple Fermi surfaces that have been smeared over momentum space? 
Strictly local criticality is known to arise in a certain limit of spin glass systems [48J. Here we 
are specifically interested in semi-local criticality, where there is a momentum dependence 
in the weights and also no entropy at zero temperature. 

We suspect that the decoupling of the gauge invariant perturbations we have performed 
may be simpler to understand by working directly in an uplifted theory with fewer basic 
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fields. These could be the uplifts in [32J or those in |27J. 

An interesting feature in both classes of holographic models that we considered is that 
there is a stronger momentum space structure in the transverse channel than the longi- 
tudinal channel. This is the opposite to free fermions. Also of interest is the fact that 
holography makes manifest the physical distinction between low energy spectral weight and 
the static susceptiblity (which controls the existence of Friedel oscillations) . We have found 
strong momentum-space features in the spectral weight. However, the static susceptibility 
either does not exhibit any features - in the probe brane case - or is independent of IR data 
- in the semi-locally critical case. We noted that the probe brane case is also known not to 
exhibit quantum oscillations as function of magnetic field. It appears to be an interesting 
question for future work to look for quantum oscillations in the semi-locally critical cases 



with < rj < 2; it is possible that the dramatic spectral weight (24) will have consequences 
for the magnetic susceptibility. The first step is to construct the extremal dyonic black 
holes in these theories. 

We end with a few comments on the D3/D5 theory. From a weakly coupled worldvolume 
perspective this theory has both charged bosons and fermions. The momentum space 
structure that we have nonetheless found in the strongly interacting gravity dual suggests 
that it is not necessary to use purely fermionic worldvolume theories to obtain strongly 
interacting cousins of Fermi liquids. Because we have treated the D3/D5 system in the 
probe limit, the Drude delta function is not visible in the conductivity [H]. It follows 
that we cannot compute the d.c. conductivity due to charged lattices or impurities in this 
case, as these revolve around resolving the Drude delta function [H] . It would certainly be 
interesting to find a setup where the momentum space structure we have uncovered, which 
have the form of smoothened out Fermi liquids, could be used as an input for the study of 
d.c. transport. 
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A Free fermions 



The propagator of a spinless nonrelativistic free fermion in Euclidean space is 

G(iu n ,k) = . 1 , (70) 

with 

&-£-*». Pi) 

Here m is the mass and \x the chemical potential (Fermi energy). The Euclidean density- 
density Green's function is, not keeping track of the overall normalization, 

(J* J*) (iw„, k) oc y dJ n d 2 k'G(iuj n + iu' n ,k + k')G(iu)' n , k') , (72) 

oc /wfhkh?!^), (73) 

oc [ F qdq [ d6 [ p^— r J .(74) 

Jo Jo \iu)n- m - ^cos9 iun+2Wl-^cos8j 

The retarded Green's function is given by taking iu n — )■ u + ie. The imaginary part is now 

particularly easy to compute 

lmG$ Jt (a J ,k)^^±j\d q £\e5^ T £ i - 1 ^co S e"j . (75) 

This integral is trivial to perform and in the low energy limit we obtain 

1 

Im J j y ' oc { 

UJ 



k < 2kF 



kyjk 2 - k 2 /4 , ( 76 ) 

k>2k F 



as quoted in the introduction. 



The transverse spectral weight is similarly obtained. The Euclidean Green's function is 

(J ± J ± )(iu n ,k) oc J doo' n d 2 k' (k' 2 - ( k k f )2 ^G(iu n + iu' n ,k + k')G(iLo' n ,k'). (77) 

Following the same manipulations as in the density channel, we obtain that the imaginary 
part is given by 

ImG?, 7| (u,k) oc V± [ F q 3 dq [ sin 2 6d06 ( u T — - — cos 6 V (78) 
^ Jo Jo \ 2m m ) 

These integrals are also easily performed, and give the result for the transverse channel that 

we quoted in the introduction. 

The static susceptibilities are computed using similar methods. Taking the real instead 



of the imaginary part of the last line (74) above, after analytic continuation ioo n —¥ uj + ie, 
gives the principal value of the pole rather than a delta function. These integrals can be 
performed to give the well known two dimensional Lindhard functions. 
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B Gauge invariant variables for semi-local models 



Even though we have chosen to work in the radial gauge (A r = g^ r = 0), there is still 
some gauge redundancy in each set of variables, arising from the U(l) gauge invariance and 
diffeomorphism invariance. Invariance under a U(l) gauge transformation parametrized by 
A and a diffeomorphism generated by the vector field £^ tells us that the modes 

= Sg^v - ^g^v = Sguv - (V M ^ + V^) , (79) 
5A^ = - £ e i, - V,A = - (f V.i, + i.V.f ) - V,A , (80) 
6$ = 5^-C^ = 5^-^V^, (81) 

represent the same physical perturbation as {5g^ u , SA^, 5&}. Here C^X is the Lie derivative 
of X along £^ and we have denoted quantities associated with the background with bars 
(i.e. V|» is the covariant derivative associated with the background metric g^v). 

With Fourier transformed fluctuations {Sg^u, SA^, 5&} and parameters {A, of the 
form X(t,r,x) = e'^-^X^r), the game is to find linear combinations of the fluctuations 
that are invariant under the transformations above while staying within the radial gauge 
and staying within the correct parity channel. In the transverse channel, the three modes 
{SA y , 5g xy , Sg y t} can be combined to form the following two gauge invariant variables: 

V>_l,i = 5 A y , ^±,2 = r 2 (uj5g xy + k5g yt ) . (82) 

In the longitudinal case a suitable set of gauge invariant modes is 

, K r 3 h'(r)d r [f(r)g(r)]-2f(r)g(r)h"(r) 1 
%i = ^rSA t - — L2/(r)g(r) S 9vv + ^f{r) h (r)5g tt (83) 

= —Q SA t - ^ l + r l § 
eL ZoL 2 r]r 1 ^ 2 ^x( r ) ' 

^ = S9 " + 2r,(r) + fc4,4)+rV(r) ar[r2(5fe + M (M) 
= §9yy + ^(l + ,)(l + 2F)-(l + 2,) X (r) 5r ^ 2 {59XX + 59vy)] ' 



r 3 <J>'(r) . 
6 + ~2i^ s to» 



5§ r 



2 



71 + 27^ 2L2 



+ 7m 5 9yy > ( 85 ) 



where we have used the background solution fel) and defined fc = yT+ 1/rjk/Vo to simplify. 
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C Perturbation equations for semi-local models 



C.l Transverse modes 



The linearized Einstein-Maxwell-dilaton equations of motion for the transverse gauge in- 
variant modes are found to be 



+ 



rr/x(r) 

4(1 + 77) 2 [r 4 /^o;2 _ k 2 X (r)] 
r 2 VQT] 2 x 2 {r) 



(86) 



8w 2 (i + n) 



r 2(l-2/r,) x ( r ) ^4/^2 _ fc2 x ( r )] 



^±,l( r ) ' 



, fc 2 (7?-2) X 2 (r)+r 4 /^ 2 [r ?X (r)-2(l+r ? )] 4fcr 1 ~ 2 /" ? , 
= ^,2^) + ; x r ,/„ 9 to / m ^±,2(0 " ^±,l( r ) ( 87 ) 



r?7x(^)[^ 4 / T 'w 2 — & 2 x( r )] 



8fcr 2 /^ 2 [(r ? -l) X (r)-(l + 7 ? )] 7 4(1 + r?) 2 ^ 4 /^ 2 - fc 2 X (r)] 

2 7 U 4/ n 9 7 2 I M V>±,l( r ) + 2t/2 2 21 \ V±M r ) > 

where we have redefined ^±,1 to include a trivial constant: l = Zq^JI + rjifrx i. Taking 
our cue from [19], these equations can be decoupled by introducing 

r 1/2+1/57^3/2^) 



where 



V-L+0) 



2r l-2/77 



r 4 / 7 ^ 2 — fc 2 x( r ) L J 



PiM^-L.lO*) +^±,2(0 



(88) 



p±(r) 



( r 4/r? w 2 _ k 2 x (r)) ^-1 + a/i + 27?fc 2 ^ - 2fc 2 x( 



(89) 



fe»7X(r) 

The equations of motion for these new variables are then the decoupled second order equa- 



tions 



= <p" (r) I ( q±[xir)] I 4 ( 1 + ^) 2 [^ 2 -fc 2 x(0] V^ ±(r) 

t, I ^^2 r 2^2^ r ) ^A2^2 r 2^,2^ r ) I 



(90) 



where <?±[x( r )] is given by: 



q± [x(r)\ 



4(1 + r\Y - 8(1 + 77) 3 - ry ± 2v/l + 2r?A; 2 x(0 - 15^x W 



(91) 



C.2 Longitudinal modes 



The linearized Einstein-Maxwell-dilaton equations of motion for the gauge-invariant longi- 
tudinal modes are found to take the form 

= i>Ur) + mArWyir) + % W| U (r) + 1 ^2^) ~ % M > (92) 
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where a 12 = «2i = «3i = «32 = 012 = 032 = 0; 
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0, 



an 



«22 



«23 



-2(1 + 77) + (4 + 7r])x(r) _ 1 _ 2(1+7/) 

rrjx(r) r rrj X (r) 

r ~ rl X \t)t) { {l + ^ 2(1 + 2P) " (1 + ^ )(4 + 5??)(1 + 2PM? 

+ [3 + 5r/ + 4r/ 2 + 4(1 + r/)(2 + t/)P] x 2 (r)} , 
8(l + T7) 2 (2 + T7)(l + 2fc 2 )( X (r)-l) 



and 

0n = 

021 



r 3 r/7(r) 

-2(1 + 77)[12 + r/(16 + 377)] + (2 + t?)(2 + 3r/) 2 X (r) 



r 2 77 2 (2 + r/)x(r) 



013 = 



. 16 p r -5-2/^( 1 + r? )2 



r? 2 x ( r ) 



-4(1 + 7/) 



{(1 + 7/)(4P - 1) - 2t/(1 + 2P) X (r) + 3(7/ - l)x 2 (r)} , 



r i 2 / r >r]x(f)'y(,f) 

i 22 = — ^ - — (2(1 + 7/) 3 (2fc 4 - P - 1) - 2(1 + 7/) 2 (l + 2fc 2 ) [(2 + r/)P-3(l + 77) X (r) 



r2r //C( r )7( r ) 



-(1 + 7/) 4 + 5r/ + 6t/ 2 + 2(5 + 67/ + 4r/ 2 )P xV) 
+ 7/(r/-l)(l + 2r/)x 3 (r)}, 



23 



-8(1 + 7/) 







31 



r 4 r?x(r)7(r) {(! + ^) 2 ( 4p - !) - 2(1 + r7)(l + 2P)[(2 + t/)P - l]x(r 

+ (7/-l)[l + r/ + 2(2 + 7/)P]x 2 (r)}, 
-4r 1+2 / , '(l + 7/) 



7/(2 + r))x(r) 
with fc = y/l + I/77 fe/Vo and 7(7*) given by 

7 (r) = [(1 + r/)(l + 2fe 2 ) - (77 - l) X (r)l [(1 + r?)(l + 2k 2 ) - (1 + 27/) X (r) 



(93) 



To clean this up, we first get rid of the first derivative terms by writing the modes as 
= Ci(r)^\\,i, where the Q(r) are given by 

£ ( r ) - r ~ 5/2 ~ lA; c (r) _ r- 3 / 2+lA ? (1 + r/)(l + 2fc 2 ) - (7/ - l)x(r) ^ r V2+iA; 



y/xTO (l + 7/)(l + 2p)-(l + 2r/)x(r-) 



With these definitions, the modes ip\\,i/3 satisfy equations of the form = ip'^ + Pij^j, 
with = 032 = 0. The equation for ipu 2 is n °t of this form, but if we look at the linear 
combination ip^ 2 = ip^ 2 + J'( r )V'||,3> where 

X(r)-1 (l + r/)(2 + 7/)(l + 2fc 2 ) 



p(r) 



1 + (1 + ri)k 2 (1 + t/)(1 + 2P) - (7/ - l)x(r) 



(94) 
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then <£>|| j2 satisfies an equation of the form = ip" i + fajip^j. Note that at T = 0, x( r ) = 1 
so p(r) = and this combination is unneccesary to decouple the fields (e.g. a 2 3 = at 
T = 0). Finally, we clean up notation a bit by letting ip^ = tp^^ for i = 1,3 with c£>|| ;2 
defined as above, so the equations of motion for the ipn i are 



4(1 + t/) 2 [r 4 A'w 2 - k 2 X (r) 



-<P\\,i( r ) ' 



r 2 77 2 VJ, 2 x 2 (r) 

where /3i 2 = /?3 2 = - this means that there are self-contained equations for and p\\^ 
that do not involve <^ii 2 with 



3 1 2(l+r ? ) 2 (10 + r ? ) (1 + r/) 2 

Pll — ~T^T n x / \ H 



21.2 



4r 2 r 2 r/ 2 (2 + f])x( r ) r 2 r] 2 x 2 (r) 
-4(1 + 77) 



r 2 ?7(2 + ?7)x(r) 
while ip» 2 couples to both tpn-y and </?ii 3: 
-4(1 + ??) 



013 = 



Ass = 



-16(1 + 77)^ 
r 2 r] 2 x(r) 



2(l+r?) 2 + (I + 77) 2 



#21 = 



r 2 7?x(r)7(r) 



'22 



|(1 + f])(l + 2fc 2 ) 2(1 + r])k 4 - (1 + 3r?)P - 2 - 7/ 
+ 2(1 + 2P) r ? 2 + (l + r ? )P x(r) + (r? - 1) 2 - 77 + (1 + # * 2 (r)} > 

4^x4)^) ( 4(1 + 77)4(1 + 2P)2 + 16(1 + ^ + i 2 " ^ + (3r? " 1)P ] XW 

-(I + 77) 2 76 + 160fc 2 + 8r/(l + 4P) 2 + r/ 2 (124fc 4 + 28P-21) xV) 

- 2(r/ 2 - 1) [l2 + I677 + 5r/ 2 + 2 (16 + 77(16 + 77)) fc 2 l x 3 (0 
+ 77 2 ( 77 -l)V(r)}, 

&3 = ^(^ {(1 + V) 2 (l + 2P) [277(1 + # - 2(77 2 - 1)P - 77; 

- 2(1 + 7?)(1 + 2k 2 ) [77(1 + 7?)(2 + 7?)P 

+ (2(l + 7/)- 77 2 (2 + 77)) P-77] X (r) 

- (7/ - 1) [2(1 + t?)(2 + 7/)fc 4 + 2k 2 - 7/(7/ + 1)] X 2 (r)} • 



Here 7(f) is given by 



7(r)= l + (l + 77)fc 2 (l + 77)(l + 2fc 2 ) - (r/- l)x(r) 



(95) 



The T = case simplifies considerably, and we can fully decouple these equations; we do 
this in the main text. For T 7^ 0, <p^i and </?|| i3 are easily decoupled. In general, it is not 
clear whether a decoupled equation can be found for p\\^, except in the case 77 = 1 for which 
simplifications occur, as we discuss in the main text. Even if we were able to decouple the 
equations for T / and 77 / 1, we don't expect to be able to solve them analytically. 
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